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A  METHOD  TO  COMPUTE  THE  CONTACT  FORCE  OF  A  BODY 
IMPACTING  A  RING-STIFFENED  CYLINDRICAL  SHELL  USING 
A  LUMPED  PARAMETER  FINITE-DIFFERENCE  MODEL 


INTRODUCTION 

The  Naval  Research  Laboratory  has  in  the  past  beer  involved 
in  defining  shock  design  inputs  for  equipment  aboard  submarines. 
The  Navy  is  currently  involved  in  s e a -e xe r c i s e s  in  which  inert 
torpedoes  with  high  velocities  are  being  used.  With  the 

possibility  of  a  torpedo  impacting  a  submarine  and  damaging  vital 
equipment,  resulting  in  loss  of  life  and  the  submarine,  a 
reliable  method  to  estimate  the  impact  force  is  needed.  This 

report  presents  a  method  to  compute  the  contact  force  of  a  body 
such  as  an  inert  torpedo  impacting  a  r ing -s t i f f ene d  cylindrical 
shell  using  a  lumoed  parameter  finite-difference  computer 
program . 

The  method  is  very  simple  in  concept.  The  measured  velocity 
of  the  struck  region  is  impressed  on  the  lumped  mass  around  the 
impact  point  in  the  computer  model.  The  computer  program  for 

each  time  step  calculates  the  velocity  of  the  struck  mass 

assuming  no  impact  force.  All  the  internal  forces  on  the  mass 
are  calculated  by  the  program.  The  calculated  velocity  is  in 
error.  The  force  necessary  to  move  the  mass  at  the  correct 

velocity  (impressed  velocity)  is  computed  at  each  time  step. 
This  additional  force  is  the  impact  force. 

The  method  has  several  important  advantages.  The  data 
measured  at  the  impact  point  are  used  to  drive  the  mathematical 

model.  The  response  of  the  model  at  each  time  step  gives  the 
impact  force  at  the  previous  time  step,  and  for  shock  type 
loadings  of  10-30  ms  the  model  has  little  time  to  degrade  in 
accuracy.  In  addition,  because  of  the  short  duration  of  the 
calculation,  damping  plays  little  part. 

A  previous  report  "A  Method  to  Compute  the  Force  Signature 
of  a  Body  Impacting  on  a  Linear  Elastic  Structure  Using  Fourier 
Analysis”  [1]  discusses  a  frequency  domain  technique  for  deriving 
the  impact  force.  This  method  suffers  from  several  major 
drawbacks.  In  order  to  obtain  convergence  of  the  Fourier 

integrals  of  the  impulse  response  and  response  to  force, 
functions  out  to  long  times  must  be  used.  The  long  times 
required  introduce  errors  in  the  derived  force  due  to  errors  in 

the  measured  gage  response  and  errors  in  the  computer  impulse 

response,  especially  at  late  times.  In  order  to  use  the  computed 
impulse  response  damping  must  be  applied  to  the  function. 
Obtaining  the  correct  damping  as  function  of  frequency  is  another 
difficulty,  and  introduces  errors  in  the  derived  force. 
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EXPERT MEN 


TORPEDO  IMPACT 


In  1932,  a  series  of  rocket  assisted  air  drops  of  various 
torpedoes  was  conducted  by  NSRDC-UERD  ( Naval  Ship  Research  and 
Development  Center  -  Underwater  Explosion  Research  Detachment)  at 
Portsmouth,  Virginia.  The  torpedoes  were  dropped  aeainst  a  large 
ring-stiffened  cylindrical  shell  made  cf  HY-130  steel,  also  known 
as  the  A  3-1.  Data  were  measured  on  the  A  3-1  by  velocity  meters 
anc  accelerometers  during  the  impact.  The  locations  cf  these 
velocity  meters  and  accelerometers  in  the  vicinity  of  the 
impacted  frame  are  shown  in  Figure  1. 


One  of  the  itr.  pactors  dropped  was  a  simulator  torpedo  (an 
impactor  designed  to  simulate  the  force  of  an  actual  torpedo). 
Data  was  sent  to  NRL  on  digital  tape.  The  method  as  described  in 
this  report  is  applied  to  this  data  set  as  a  test  case.  The  data 
from  the  impact  region  are  used  to  generate  the  force.  This 
impact  force  is  then  used  with  the  damped  computed  impulse 
response  to  generate  predicted  responses  for  various  gage 
locations.  From  the  predicted  responses  predicted  shock  spectra 
may  be  obtained  and  compared  with  measured  shock  snsctra. 
Predicting  shock  spectra  is  of  prime  importance  in  desist  and 
analyses  of  internal  equipments  by  current  methods. 


METHOD  USING  VGSM  TO  FIND  IMPACT  FORC! 


VGSM  (the  Variable  Geometry  Submarine  Model)  [2]  employs  a 
lumped  parameter  finite  difference  approach  to  model  the  response 
of  ring-stiffened  cylindrical  shells.  The  program  solves  the 
partial  differential  equations  of  Timoshenko  beam  theory.  A 
model  generator  program  is  supplied  to  generate  the  cylindrical 
shell  with  stiffeners.  The  user  must  write  a  main  program  which 
contains  the  additional  equipment  equations  and  calls  the  dynamic 
subroutine  which  updates  (integrates  in  time)  the  response  of  the 
masses  modeling  the  shell.  The  main  program  integrates  the 
equipment  equations.  The  model  time  histories  to  be  printed  out, 
or  saved  on  disk/tape,  must  be  defined  in  the  user-supplied  main 
program . 
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where  M  is  the  diagonal  lumped  mass  matrix,  K  is  the  symmetric 

stiffness  matrix  of  the  equipment,  X  is  the  column  vector 
equipment  mass  displacements,  Y  is  the  column  vector  of  VGSM  ma 

displacements  at  the  equinmor.  t  supports  and  T  is  a  matr 

transforming  VGSM  displacements  at  supports  to  equinre 


r.  1 
o  f 
s  s 
i  x 

L  I  v- 


9 


•  vVv 


displacements  at  lunped  masses.  The  reaction  forces  at  the 
equipment  supports  can  be  expressed  as: 

M. 


i 


[FR]  =  [  T]  [XI  IX  - 


(2) 


The  measured  velocity  of  the  shell  at  the  impact  point  is 

low  pass  filtered  tc  remove  noise  and  obtain  a  smooth  signal. 
The  VGSM  mass  corresponding  to  the  impact  point  is  moved  with  the 
filtered  measured  velocity  VM_.  at  time  t^.  The  A3-1  model  and 

equipment  is  then  allowed  to*  move  for  one  time  step  and  VGSM 

returns  for  the  struck  mass  a  velocity  of  VY^  +  j  at  time  tJ+j  , 

which  differs  from  VM.+j  ,  the  measured  velocity  at  t  +  ^  .  The 

impact  force  TF  needed  to  make  the  VGSM  mass  velocity  equal  to 

the  measured  velocity  VM^+j  can  be  computed  from: 

T?  =  FR  +  nr  (  VY .  -  VM.  )  (3) 

l  l  At  l+l  i+l 


This  equation  may  be  turned  around  to  read: 

FRi-TFi 


VM  =  VY  + 
i+l  i+l 


M 


A  t 


(4) 


This  equation  says  that  the  measured  velocity  at  t^j  is  equal  to 

the  VGSM  velocity  at  t . +  ^  corrected  by  the  velocity  change  due  to 

F  R  i 

the  upward  equipment  reaction  force  ■■■  At  at  ti  and  the  downward 
•  -  TF  i  M 

impact  force  -  At  at  t..  The  VGSM  mass  velocity  is  then  set 

to  the  measured  ^velocity  at  t^  +  ^>  the  model  is  allowed  to  move 
for  one  time  step,  and  equation  3  is  again  used  to  find  TF.+j. 
The  sequence  for  finding  TF  is  programmed  into  the  user  supplied 
main  program. 


RESULTS  OF  METHOD 

The  above  method  is  applied  to  compute  the  contact  force  of 
a  simulator  torpedo  dropped  on  the  AB-1.  Figure  2  shows  the 
unfiltered  velocity  as  measured  by  VI 1,  on  the  end  of  the  flange 
under  the  impact  on  frame  9.  Figure  3  shows  the  velocity  record 
from  VI 1  low -pass  filtered  below  500Hz.  The  filtered  velocity 
V 1 1  drives  the  VGSM  struck  mass.  Figure  U  shows  the  impact  force 
calculated  at  each  time  step.  The  negative  values  represent 
modeling  errors  since  there  is  no  mechanism  for  the  imoactor  to 
pull  up  on  the  AB-1,  creating  negative  or  tensile  forces.  Figure 
5  shows  the  impact  force  with  negative  values  removed  and 
terminated  at  a  time  corresponding  to  no  additional  impulse  being 
delivered.  Figure  6  shows  the  momentum  or  impulse  delivered  to 
the  A3-1.  After  the  initial  momentum  of  4663  Ib-s  is  expended  a 
substantial  additional  iooulse  is  calculated  due  tc  a  large 
bounce  of  the  simulator  also  observed  during  the  drop.  Some 
impactors  do  not  bounce  and  the  calculated  impulse  expanded  then 
equals  MV  ,  the  initial  momentum  of  the  impactor.  Figure  7  shows 
the  velocity  of  the  simulator  center  of  gravity.  The  simulator 
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slows  down  and  Chen  reverses  direction  as  it  bounces.  Figure  3 
shows  the  displacement  or  crushing  of  the  inpactor.  The 
simulator  crushes  about  5  inches  during  impact,  close  co  the 
measured  6  inches.  Figure  9  shows  the  enersrv  expended  h'T  the 
simulator  during  impact.  The  initial  kinetic  energy  (1/2  MV  -  ) 
of  159,017  lb -ft  of  the  simulator  is  totallv  extended. 

COMPARISON  OF  PREDICTED  AND  MEASURED  RESPONSES 

One  method  of  comparison  of  ore dieted  versus  measured  is  to 
move  the  VG3M  struck  mass  with  the  measured  velocity  and  compare 
the  velocity  computed  to  the  measured  velocity  at  the  neighboring 
frames  and  on  the  simulated  equipment.  This  amounts  to  comparing 
measured  and  comouted  velocities  during  the  time  that  the  force 
is  computed.  Figures  10,  11  and  12  show  good  agreement  between 
predicted  and  measured  time  histories  on  frames  in  the  vicinity 
of  the  impact  frame.  Figures  13  and  1 R  also  show  good  agreement 
between  predicted  and  measured  velocity  time  histories  on  the 
simulated  equipment.  It  should  be  noted  that  the  measured  data 
has  an  arbitrary  zero  time,  however  by  shifting  the  measure  curve 
right  or  left  a  better  match  can  be  obtained. 

Another  method  of  comparison  is  to  convolve  the  truncated 
force  (no  negative  values  and  terminated  after  impulse  is 
expended)  and  a  YGSM  computed  impulse  response  to  obtain  a 
oredicted  response.  The  computed  impulse  response  should  be 
damped.  The  following  method  is  employed.  Shock  spectra  of 
experimental  records  tor  various  lengths  of  time  are  produced. 
The  growth  of  peaks  at  different  frequencies  in  che  shock 
spectrum  is  used  to  estimate  an  equivalent  viscous  damping  factor 
in  five  frequency  bands  from  0  to  500  Hz.  The  estimate  is  based 
on  using  the  growth  of  the  analytical  solution  of  an  undamped 
shock  spectrum  oscillator  to  damped  base  motion.  The  ratio  of 
the  amplitudes  at  two  different  times  gives  the  damping  in  the 
base.  Damping  factors  obtained  by  this  technique  ranged  from  IN 
to  3”  of  critical. 

Impulsive  responses  for  the  instrumented  frames  and 
locations  on  the  simulated  equipmer t  are  found  using  VGSM.  Each 
impulse  response  is  broken  into  the  same  five  frequency  bands  as 
the  data  by  using  band -pass  filters,  and  the  time  histories  for 
each  band  are  modified  by  the  exponential  damping  factors.  The 
damped  time  histories  for  each  of  the  five  bands  are  summed  to 
form  a  damped  impulse  response.  The  response  of  a  linear  elastic 
system  such  as  the  A  B  - 1  may  be  computed  from  the  convolution 
integral  : 

t 

R  i  (  t  )  -  /  F  (  T  )  h  i  (  t  - 1 )  d  t  (  5  ) 
o 

This  equation  is  the  convolution  of  the  impact  force  and  the 
damped  impulse  response  at  location  i.  The  R_.  (t)  are  computed 
for  locations  or.  the  impact  frame,  neighboring  f*r  a  m  e  s  ,  and  on  the 
simulated  equioment  using  a  Fourier  transform  computer  code. 
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A  valid  method  for  obtaining  the  contact  f< 
such  as  a  torpedo  impacting  a  ring-stiffened  cvlinder  has  been 
described.  The  method  is  tested  on  data  from  a  simulator 
torpedo.  The  calculated  force  is  consistent  with  measured  data 
in  many  respects;  momentum,  energy  and  crush.  The  derived  force 
in  turn  predicts  velocity  responses  and  shock  soectra  which  agree 
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